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ABSTRACT — Chromosomes and 15 different protein systems of seven taxa of the Japanese microtine 
voles were examined to establish their phylogenetic relationships, making use of differential staining and 
protein electrophoretic techniques. While interspecific variations in the size of C-bands were observed 
only in the Y chromosomes, a highly homologous G-banding pattern as well as the karyotypic similarity 
were found among the microtine taxa examined, indicating that all of the micro tine taxa here dealt with 
are regarded to be closely related to each other. The biochemical dendrogram for the seven taxa 
constructed from the Nei’s genetic distances between taxa by using the UPGMA clustering method 
showed that the seven microtine taxa could be classified into three groups, though the differences of the 
D values were small in general, and this classification was well consistent with the karyological evidence 
of the interspecific Y-chromosome variations and their geographic distribution pattern. Furthermore, 
the electrophoretic results indicated that the andersoni complex should be classified not as a member of 
the genus Aschizomys or Eothenomys, but as that of the genus Clethrionomys, and that both andersoni 
and niigatae may still be at the subspecies level from their high genetic similarity. Phylogenetic 
relationships among the Japanese microtine voles are discussed in some detail from karyological and 
biochemical viewpoints. 


INTRODUCTION 

Ever since the first description of Evotomys (= 
Clethrionomys ) smithii, E. bedfordiae, E. anderso- 
ni and E. mikado [1,2], taxonomy of the Japanese 
microtine voles has been studied by many workers 
mainly from a morphological standpoint [3-9]. 
However, in spite of the extensive examination on 
the phylogeny of these microtine voles, their phy- 
logenetic relationships are still in a great con- 
troversy. 

While the voles of Hokkaido and its adjacent 
islands have been classified into either two species 
[9-11], three species [12-14] or five species [15, 
16], the voles of Honshu, Shikoku and Kyushu into 
either two species [9, 14, 17] or five species [5, 6, 
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15, 18] excluding Microtus montebelli. On the 
other hand, four different views on the generic 
allocation of andersoni have been proposed on the 
basis of the traditional taxonomic criteria such as 
dental and skull systems: andersoni should be 
included in either the genus Clethrionomys [14, 
17], the genus Aschizomys [5], the subgenus Aschi- 
zomys of the genus Clethrionomys [6], or the genus 
Eothenomys [9]. The specific (and even generic) 
allocation of andersoni has still remained indistinct 
and unresolved. 

In the last decade, chromosome banding and 
protein electrophoretic techniques have provided 
much relevant, in some cases critical, information 
on the phylogenetic relationships in various groups 
of mammals [19-22]. As to the Japanese microtine 
species, no systematic investigation on the 
chromosome banding pattern as well as on the 
electrophoretic analysis of proteins has been made 
so far, except for a couple of case reports on the 
chromosome banding analysis [23-25]. In this 
study, we have attempted to establish the phy- 


410 


I. Yoshida Y. Obara and N. Matsuoka 


logenetic relationships among the Japanese micro- 
tine species by using chromosome banding and 
protein electrophoretic techniques. 

MATERIALS AND METHODS 

Animals 

A total of 72 specimens from seven microtine 
taxa were examined in this study. They were all 
trapped alive at various localities of Honshu, Hok- 
kaido and Rishiri Island. For nomenclature of the 
voles studied here the taxonomic system of Corbet 
[14] was adopted. Triplet and doublet of alphabets 
in parentheses show abbreviations for their scien- 
tific names. 

Clethrionomys rutilus mikado (Crm): Mt. Pete- 
gari (1050 m), Hokkaido ($0, £2); C. rex (Cre): 
Rishiri Shrine and Kanrosen, Rishiri Is. ( $4, £4); 
C. rufocanus sikotanensis (Crs): Rishiri Shrine, 
Rishiri Is. ($2, £0); C. rufocanus bedfordiae 
(Crb): Mt. Petegari (1050 m), Shojiyama, Yunotai 
and Miyauta, Hokkaido ( $ 13, £ 11); C. andersoni 
andersoni (Caa): Mt. Iwaki, Zatoishi and Ainai, 
Aomori Pref., Nagabashiri, Akita Pref., Asahi- 
kousen, Yamagata Pref. ( $ 16, -£9); C. andersoni 
niigatae (Can): Jumonjitoge (2050 m), Saitama 
Pref. ($2, -£2); Eothenomys smithi (Es): Nanza- 
wakousen and Nezame, Nagano Pref., Eiheiji and 
Heisenji, Fukui Pref. ($5, £2). All of the 
animals examined are preserved as standard 
museum specimens (skin with skull) or formalin- 
fixed specimens in the Department of Biology, 
Faculty of Science, Hirosaki University, Hirosaki. 

Chromosone preparation 

Metaphase chromosomes were obtained from 
the femoral bone marrow cells of the chloroform- 
anesthetized specimens. The procedure for 
chromosome preparation was almost the same as 
that previously described [26, 27]. For G- and 
C-band staining, the ASG [28] and BSG [29] 
methods were adopted. Analytical data on the 
conventionally-stained and G- and C-banded 
karyotypes are summarized in Table 1. 

Electrophoresis 

Five different tissues (kidney, liver, intestine, 


pancreas and skeletal muscle) which had been 
stored at — 80° C were used for the electrophore- 
tic study. The procedures for tissue preparation 
and polyacrylamide gel electrophoresis were 
almost the same as those described by Matsuoka 
[30, 31]. The following 14 different enzymes and 
general proteins (non-enzymatic proteins) were 
assayed, using supernatants of tissue homogen- 
ates: a-glycerophosphate dehydrogenase (a- 

GPDH), glucose-6-phosphate dehydrogenase 
(G6PD), hexose-6-phosphate dehydrogenase 
(H6PD), lactate dehydrogenase (LDH), malate 
dehydrogenase (MDH), octanol dehydrogenase 
(ODH), 6-phosphogluconate dehydrogenase (6- 
PGD), xanthine dehydrogenase (XDH), superox- 
ide dismutase (SOD), hexokinase (HK), aspartate 
aminotransferase (A AT), esterase (EST), amylase 
(AMY), leucine amino peptidase (LAP), and 
general protein (GP). G6PD, H6PD, LDH, 
ODH, 6-PGD, HK and A AT were assayed with 
the extract of the kidney; a-GPDH, MDH, XDH 
and SOD with that of the liver; AMY and LAP 
with that of the intestine and pancreas, and EST 
and GP with that of the skeletal muscle. As three 
enzymes (EST, LAP and AMY) and GP showed 
many polymorphic bands, these were assayed us- 
ing a vertical slab gel apparatus (gel plate mea- 
sured lx 160x140 mm with 8—20 sample slots). 
For other enzymes a disc gel apparatus (3 X 11 X75 
mm plastic column) was used. The stain mixtures 
for LDH and 6-PGD were prepared according to 
Shaw and Prasad [32] and that for a-GPDH 
according to Ayala et al. [33]. Stain recipes for 
other enzymes were those described previously 
[31, 34, 35], GP was stained with Coomassie 
Brilliant Blue R-250 by the method of Matsuoka et 
al. [36]. 

RESULTS 

Karyological analysis 

All of the seven taxa of microtines examined had 
the diploid number of 56 (Table 1). A convention- 
al karyotype of Caa is presented in Figure la as a 
standard karyotype representing these microtines. 
The karyotype consisted of 26 pairs of acrocentrics 
of gradually decreasing size, one pair of the smal- 
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Table 1. Analytic and cytologic data on the seven taxa of the Japanese red-backed voles examined in 
the present study 


No. of cells karyotyped Sex chromosome 

Species NFA 2n 

Conv. ASG BSG X Y 


Clethrionomys 

rutilus mikado 

8 

5 

5 

56 

56 

A 

M’ 

C. rex 

23 

16 

8 

56 

56 

A 

A 

C. rufocanus 
sikotanensis 

13 

14 

6 

56 

56 

A 

A 

C. r. bedfordiae 

48 

57 

24 

56 

56 

A 

A 

C. andersoni 
andersoni 

65 

41 

30 

56 

56 

A 

M 

C. a. niigatae 

18 

15 

14 

56 

56 

A 

M 

Eothenomys smithi 

15 

12 

11 

56 

56 

A 

St 


* Shimba et al. [41], Tsuchiya and Yoshida [44] and Tsuchiya [37]. A, acrocentric; M, metacentric; St, 
subtelocentric; Conv., conventional staining; ASG, G-banding; BSG, C-banding; FNA, fundamental 
number of autosomes; 2n, diploid number of chromosomes. 


lest metacentrics and the XY chromosomes. Small 
but easily detectable short arms were found only in 
the pair No. 3 which can be regarded as 
acrocentric judging from its arm ratio, 9.87 + 0.71, 
though Tsuchiya [37] thought this chromosome to 
be subtelocentric. Thus, the NFA of Caa should 
be 56. Similar pattern of the small short arms was 
observed in all of the remaining taxa. Their 
autosomal constitution was almost identical with 
that of Caa, and hence their NFA was 56 in all the 
taxa examined (Table 1). The interspecific differ- 
ence in the length and morphology of chromo- 
somes was found only in the Y chromosomes. The 
X chromosomes were all acrocentrics of the same 
length. 

A G-banded karyotype of Caa is presented in 
Figure lb. Chromosomes are arranged according 
to the numbering system proposed by Gamperl 
[38] who has described the G- and C-banding 
patterns of the grey red-backed vole, C. rufocanus, 
and the bank vole, C. glareolus. As clarified in the 
previous report [25] , the G-banding pattern of Caa 
is almost identical with that of these continental 
species. So, making this karyotype of G-bands 
standard, detailed pair matching analysis of G- 
banding pattern was made among the seven taxa 
(Fig. 2). As clearly demonstrated in the composite 
karyotype, almost perfect G-band homology was 
obtained in all the chromosome complements in- 
cluding the X chromosomes among the seven taxa, 


reflecting their phylogenetic kinship. The C-bands 
of Caa were all centromeric except for the Y 
chromosome which was entirely heterochromatic 
along its arms (Fig. lb, insertion). There was no 
significant difference in the G- and C-banding 
patterns among the seven taxa, except for their Y 
chromosomes. The Y chromosomes showed dis- 
tinct interspecific variations in their G- and C- 
banding patterns, but no detectable variation at 
the intraspecific or subspecies level (Fig. 3). On 
the basis of the size and morphology of the Y 
chromosomes the microtine taxa under present 
study could be classified into three groups as 
follows: andersoni- type carrying a metacentric Y 
(Caa and Can), rufocanus- type carrying an 
acrocentric one (Crb, Crs and Cre) and smithi- type 
carrying a subtelocentric one (Es). The rates of 
the Ys to the Xs (Y/X + SD) were 0.28 + 0.02 in 
Cre, 0.28 + 0.03 in Crs, 0.30 + 0.03 in Crb, 0.28 + 
0.03 in Caa, 0.27 + 0.05 in Can and 0.41 + 0.05 in 
Es. Thus, the Y chromosomes correspond to 27— 
30% of the X chromosomes in length except for 
Es, in which the Y was about 41% of its X 
chromosome length. 

Electrophoretic analysis 

From the electrophoretic patterns of 14 different 
enzymes and general protein observed in and 
among the seven taxa of microtines, we assumed 
26 genetic loci (Table 2). The major features of 
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Fig. 1. Conventional (a) and G-banded (b) karyotypes of a male Japanese red-backed vole, Clethrionomys a. 
andersoni. 


variation in these proteins are summarized as 
follows. 

Eight enzymes: a-GPDH, G6PD, H6PD, ODH, 
6-PGD, XDH, HK and AAT exhibited a single 
band of activity and they were all monomorphic. 


The bands of individual enzymes showed the same 
electrophoretic mobility in the seven taxa. 

LDH showed five-banded pattern in all taxa. As 
this enzyme is known to be a tetrameric protein, 
the five bands could be interpreted as five different 
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Fig. 2. Pair-matching of G-banded haploid chromosomes from each of seven taxa of the Japanese red-backed voles. 
From left to right; Clethrionomys rutilus mikado, C. rex, C. rufocanus sikotanensis, C. r. bedfordiae, C. andersoni 
andersoni, C. a. niigatae and Eothenomys smithi. 


tetramers consisting of two polypeptides produced 
by two codominant alleles at two different loci 
(LDH-1 and LDH-2). Each of these bands showed 
the same mobility in all taxa and the two loci were 


monomorphic in each taxon. 

MDH showed two active bands, of which the 
faster band was markedly high in activity. These 
two bands were interpreted as the products of two 
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Fig. 3. G-banded (G) and C-banded (C) sex chromosomes, X and Y, of the seven taxa of the Japanese red-backed 
voles. 


different loci (MDH-1 and MDH-2) in the light of 
other electrophoretic studies of various animal 
groups. Each of these two bands also showed the 
same mobility in all of the taxa. The two loci were 
also monomorphic in each taxon. 

SOD consisted of three bands with the highest 
activity in the fastest band. These were interpreted 
as the products of three different loci (SOD-1, 
SOD-2 and SOD-3) in the light of the elec- 
trophoretic studies of sea-urchins [31, 34, 35]. The 
bands of SOD-1 and SOD-2 in Es showed the 
faster mobilities than those of other taxa, but the 
band of SOD-3 had the same mobility in all the 
taxa. These three loci were monomorphic in each 
taxon. 

EST activity was detected as several bands 
which were grouped into two zones. The slow 
zone (EST-1) consisted of a single faint band and 
was monomorphic. While the fast zone (EST-2) 
showed high genetic variability in Crb, Crm and 
Es, it was monomorphic in the remaining four 
taxa. 

LAP showed only a single active band in Crm, 
but in the remaining six taxa the enzyme consis- 
tently appeared as two active bands. These band 
patterns were interpreted to be controlled by two 
different loci (LAP-1 and LAP-2). LAP-1 was 
polymorphic in four taxa (Crb, Cre, Caa and Es) 
and LAP-2 in three taxa (Crb, Caa and Es). 

AMY activity was detected as one to two bands. 


It was polymorphic within each of four taxa (Crb, 
Caa, Can and Es), and furthermore it also showed 
considerable variations among taxa. 

In GP, six to ten bands were obtained in each 
taxon. These were assumed to be products of six 
different loci (GP-l — GP-6) from the elec- 
trophoretic mobilities of these bands. Four of 
them (GP-1, GP-2, GP-4 and GP-6) could be 
presumed to be commonly possessed by all taxa, 
but the others (GP-3 and GP-5) were considered to 
be specific only to certain given taxa. Polymorphic 
band patterns were observed in GP-3 of Crs, in 
GP-5 of Crb and in GP-4 of Caa and Can. 

The allele frequencies for all loci in these seven 
taxa are given in Table 2. In this table, it was 
assumed that the bands showing the same elec- 
trophoretic mobility between taxa in a given locus 
are products of the same alleles. Based on these 
data, the genetic identity (I) and genetic distance 
(D) between each taxon were calculated by the 
method of Nei [39]. Table 3 represents the ma- 
trices of I and D values between all pairs of the taxa 
examined. As evident from this table, the I value 
between Caa and Can, which allopatrically inhabit 
the hilly countries of northern Honshu (Tohoku 
district) and the high mountaines (more than 2,000 
m above the sea level) of central Honshu (Kanto 
and Chubu districts), respectively, is the highest (I 
—0.951). The I values between Crb,Crs, Cre and 
Crm, all of which are distributed in Hokkaido and 
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Locus Crm Cre Crs Crb Caa Can 


ff-GPDH a a 

G6PD a a 

H6PD a a 

ODH a a 

6-PGD a a 

XDH a a 

HK a a 

AAT a a 

LDH-1 a a 

LDH-2 a a 

MDH-1 a a 

MDH-2 a a 

SOD-1 a a 

SOD-2 a a 

SOD-3 a a 

EST-1 a b 


b c 


a 


c 


a 

a 

a 

a 

a 


a 


c 


EST-2 


LAP-1 

LAP-2 

AMY 

GP-1 

GP-2 

GP-3 

GP-4 

GP-5 

GP-6 


b(0.50) e b 

e(0.50) 

f d(0.29) f 

e(0.71) 

b b 

b e h 


a(0.02) e e 

W0.38) 

c(0.08) 

d(0.06) 

e(0.44) 

f(0.02) 

a(0.05) b(0.86) b 

d(0.62) d(0.14) 

f(0.33) 

b(0.43) a(0.14) b 

c(0.57) b(0.86) 

f(0.39) a (0.47) a(0.50) 

g(O.ll) c(0.53) c(0.50) 

h(0.50) 


b 

a 

a 

a 


b 


d 

a 


a(0.50) 

b(0.50) 


b 

a 


b(0.04) 

c(0.96) 


a(0.79) 

b(0.21) 


a(0.50) 

b(0.50) 


Es 


b 

b 


c(0.711 

(0.29) 


c(0.29) 

e(0.71) 


c(0.14' 

d(0.43 

f(0.43) 

d(0.57) 

e(0.43) 

b(0.21 

d(0.29) 

e(0.50) 


d 

a 


Alleles are correspondingly lettered from “a”, this being the allele of lowest mobility. The value in parenthesis 
represents the frequency of each allele in taxon. The scientific name of each taxon is shown as abbreviation 
described in Materials and Methods. 


its neighbouring islands, are relatively high with 
the range of 0.803—0.895. On the other hand, the 
I values between Es and the six taxa of the genus 
Clethrionomys are in the range of 0.713—0.792 
and lower than the values mentioned above. Fig- 
ure 4 shows the biochemical dendrogram for the 


seven taxa constructed from the genetic distance 
matrix of Table 3 using the unweighted pair-group 
arithmetic average (UPGMA) clustering method 
of Sneath and Sokal [40]. The biochemical den- 
drogram indicates the following: (1) The seven 
taxa studied here are divided into three clusters. 
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Table 3. Genetic identities (above diagonal) and genetic distances (below diagonal) between the seven 
taxa of microtines 


Taxa 

Crm 

Cre 

Crs 

Crb 

Caa 

Can 

Es 

Crm 

— 

0.803 

0.824 

0.844 

0.769 

0.713 

0.736 

Cre 

0.219 

— 

0.842 

0.855 

0.823 

0.771 

0.789 

Crs 

0.194 

0.172 

— 

0.895 

0.816 

0.763 

0.713 

Crb 

0.170 

0.157 

0.111 

— 

0.891 

0.829 

0.783 

Caa 

0.263 

0.195 

0.203 

0.115 

— 

0.951 

0.792 

Can 

0.338 

0.260 

0.270 

0.188 

0.050 

— 

0.729 

ES 

0.307 

0.237 

0.338 

0.245 

0.233 

0.316 

- 


The scientific name of each taxon is shown as abbreviation described in Materials and Methods. 


DIVERGENCE TIME (T) 
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0 (xIO 6 ) 
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Can 

Caa 

Crb 
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Cre 
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Es 


GENETIC DISTANCE (D) 

Fig. 4. A biochemical dendrogram showing the phylogenetic relationships among the seven taxa of the Japanese 
microtine voles. 


The first cluster consisted of two andersoni subspe- 
cies (Caa and Can), the second cluster four taxa 
(Crb, Crs, Cre and Crm), and the third cluster only 
Es. (2) Of the seven taxa, Caa and Can of the first 
cluster is the most closely related to each other (D 
=0.050). (3) In the second cluster, Crb and Crs 
are the most closely related to each other (D = 
0.111), and Cre is more closely related to Crb and 


Crs than Crm. However, the differences of the D 
values are small in all pair-groups of these four 
taxa. (4) Es is slightly differentiated from the first 
and second clusters, and the mean D value be- 
tween Es and two large clusters consisting of six 
taxa of the genus Clethrionomys is 0.279. 

The proportion of polymorphic loci (P) in each 
taxon was as follows: 4.0% for Crm and Crs, 3.8% 


Phylogeny of Microtine Voles of Japan 


417 


for Cre, 19.2% for Crb, 16.7% for Caa, 8.3% for 
Can and 20.0% for Es. The mean value of P was 
10.9%. The expected mean heterozygosity per 
locus (H) in each taxon was as follows: 2.0% for 
Crm and Crs, 1.6% for Cre, 8.9% for Crb, 5.5% 
for Caa, 4.2% for Can and 10.2% for Es. The 
mean value of H was 4.9%. 

DISCUSSION 

The karyotypes of the seven taxa of microtines 
here studied were essentially identical, except for 
the Y chromosome of Can, to those of the previous 
observation [37, 41-44]. The Y chromosomes of 
the andersoni complex were metacentric in both 
Caa abd Can, as far as the present specimens are 
concerned. Hsu and Benirschke [42] regarded the 

Y chromosome of Can as the smallest metacentric 
element, but Tsuchiya [37] as an acrocentric one. 
Our finding is, therefore, consistent with the for- 
mer view. These fects may be suggestive of the 
polymorphic nature of the Y chromosome of this 
taxon, just indicating an additional case of the 
parallelism of the Y chromosome variabillity re- 
ported in three Clethrionomys species; C. 
glareolus, C. rutilus and C. rufocanus [45]. 
However, it is still uncertain whether the Y 
chromosome variants have been established as the 
interpopulation polymorphism or not. Further 
examination of the specimens from various locali- 
ties is necessary for elucidation of this subject. The 

Y chromosome of Cre showed a similar but some- 
what indistinct variation: it was regarded as a 
subtelocentric element by Tsuchiya [37], and as an 
acrocentric one in the present specimens. Further, 
he is of opinion that the X chromosomes of Cre, 
Caa and Es are subtelocentric. But they were all 
acrocentric in the present study (Fig. 3). These 
minor discrepancies may be not essential but at- 
tributable to the technical matters or how to set up 
the criteria for chromosome morphology. 

As shown by the composite karyotype of Figure 
2, the seven taxa studied are almost identical to 
each other in their G- and C-banding pattern as 
well as in their conventional karyotypes, though 
their Y chromosomes showed three types of inters- 
pecific variations; andersoni- type, rufocanus-type 
and smithi-type (Table 1 and Fig. 3). These 


findings plainly signify that all of these microtines 
are the descendants from a common ancestral 
form, and with the only exception of Y chromo- 
some, no karyotypic differentiation has proceeded 
in the course of phenotypic differentiation. As 
clearly demonstrated from the pair-matching 
analysis of G-bands between Caa and C. rufocanus 
from near Gallivare, Sweden [38], the red-backed 
voles of Japan could be closely related to the grey 
red-backed vole, C. rufocanus which is widely but 
sporadically distributed from Scandinavia to Siber- 
ia, Sakhalin and Hokkaido. In view of these facts, 
it may be most likely that the seven taxa of 
microtines studied have derived from a rufocanus- 
like ancestor of the Eurasian continent. 

The biochemical dendrogram (Fig. 4) demon- 
strated that the seven taxa of microtine voles are 
genetically divided into three groups. This rela- 
tionship is well consistent with their distribution in 
Japan: Caa and Can are distributed in northern 
and central Honshu, Crs, Crb, Cre and Crm are 
endemic to Hokkaido and its neighbouring islands, 
and Es inhabits Kyushu, Shikoku and southwest- 
ern Honshu. Further, the biochemical dendro- 
gram consisting of three clusters is also consistent 
with the morphological variations of their Y 
chromosomes: metacentric in the first cluster, 
acrocentric in the second one and subtelocentric in 
the third one. The present electrophoretic results 
seem to be compatible, in a broad sense, with the 
taxonomic system of Imaizumi [15] based on the 
phenotypic characters. He classified, in the light of 
the rooting pattern of molars, these microtine taxa 
into three different genera: the genus Aschizomys 
including Caa and Can, the genus Clethrionomys 
including Crs, Crb, Cre and Crm, and the genus 
Eothenomys including Es. However, this clas- 
sification is out of harmony, in a strict sense, with 
our taxonomic system in some respects. He re- 
garded andersoni and niigatae as two distinct spe- 
cies of the genus Aschizomys. As already men- 
tioned, the D value between the two taxa is 
markedly low (D =0.050) and comparable to the D 
values observed between conspecific geographic 
populations in many other animals [46-48]. Judg- 
ing from their highly close similarity at molecular 
and chromosomal level, andersoni and niigatae 
should be considered as two subspecies of the same 
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species. Further, the mean D value (D — 0.229) 
between the andersoni complex (Caa and Can) and 
the four taxa of red-backed voles (Crs, Crb, Cre 
and Crm) is much lower than the D values re- 
ported between different genera of many other 
animals and rather at the subspecies or closely 
related species level [46-48]. Since they are also 
very similar to each other in phenotypic and 
karyological aspects, these micro tine taxa should 
be considered as the congeneric members of 
Clethrionomys, as proposed by Corbet [14]. 

In the cluster consisting of the rufocanus com- 
plex (Crb and Crs), Cre and Crm, Crb has the 
closest affinity with Crs and the D value between 
them is the lowest (D = 0.111), which is compara- 
ble to the D values observed between subspecies of 
many other animals [46-48]. Taking their distribu- 
tion areas and high genetic similarity into consid- 
eration, it seems very reasonable to regard the 
insular (Rishiri Island) taxon Crs and the mainland 
(Hokkaido) taxon Crb as conspecific, though these 
two have been regarded as distinct species by 
Imaizumi [15]. It is clear from the electrophoretic 
results that Crm had already differentiated prior to 
the evolutionary divergence of the rufocanus com- 
plex and Cre. These biochemical findings may 
substantiate the phenotypic evidence that Crm 
significantly differs in the molar pattern as well as 
in the body size from rufocanus complex and Cre. 
However, the high I values found between these 
four taxa strongly suggest that rufocanus, rex and 
rutilus are distinct but closely related species each 
other. 

With respect to the taxonomic allocation of the 
andersoni complex (Caa and Can), Aimi [9] pro- 
posed that andersoni (including niigatae) should be 
referred to the genus Eothenomys on the basis of 
the craniometry and the occlusal pattern of molars. 
If andersoni and smithi are congeneric with each 
other, then rufocanus, rex and rutilus also must 
necessarily be included in the genus Eothenomys, 
judging from the biochemical dendrogram shown 
in Figure 4. But,the view that the rufocanus- group 
must be included in the genus Eothenomys may be 
hard to accept even from usual taxonomic criteria 
because of the clear difference in rooting of the 
cheek-teeth. On the whole, our classification 
system as to the generic allocation of the andersoni 


complex is consistent not with that of Aimi [9], but 
with that of Corbet [14] who classified the ander- 
soni complex as a member of the genus 
Clethrionomys. According to Corbet [14], the 
Eothenomys species closely resemble 
Clethrionomys, only being distinguished mainly by 
the absence of rooting of the cheek-teeth even in 
old age. Further, the Eothenomys species are 
distributed mainly in China, and quite rarely in the 
Korean Penninsula and Japan, and the 
Clethrionomys species are widely distributed to the 
northern regions of the Eurasian Continent. Thus, 
it would be most probable that the genus Eothe- 
nomys might have diverged from the red-backed 
voles, Clethrionomys, by modifying the rooting 
pattern of molars somewhere in the Eurasian 
Continent and the ancestral population of Es 
might have migrated into Japan through a south- 
ern route, or the Korean Penninsula. The genetic 
differentiation of the genus Eothenomys from the 
lineage of the genus Clethrionomys may well be 
reflected in our biochemical dendrogram. 

In general, biochemical dendrogram shows not 
only the phylogenetic relationships, but also the 
sequence of evolutionary divergence. According 
to Nei [49], genetic distance (D) corresponds well 
with the divergence time (T) from the common 
ancestor, and T of two taxa can be estimated by T 
=5xl0 6 D (year). Applying this equation to our 
biochemical dendrogram, each divergence time 
may be calculated as follows: 1.4 million years 
(MY) for Eothenomys and Clethrionomys, 1.1 — 
1.2 MY for the two large groups of Clethrionomys, 
0.5 — 1.0 MY for the four taxa of Clethrionomys 
endemic to Hokkaido, and 0.25 MY for the two 
andersoni subspecies. 

For further clarification of the phylogenetic rela- 
tionships among the Japanese microtine voles, it 
would be needed to synthetically examine the 
chromosomes, genetic distances and the phenoty- 
pic characters of all the taxa including C. monta- 
nus, C. imaizumii and E. kageus which could not 
be dealt with in this study. 
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